The Statistical Energy-Frequency Spectrum of 
Random Disturbances 

By JOHN R. CARSON 

A mathematical discussion of the statistical characteristics of Random 
Disturbances in terms of their " energy-frequency spectra " with applica- 
tions to such typical disturbances as telegraph signals and " static ". 

IN a paper entitled "Selective Circuits and Static Interference" 
(B. S. T. J., April, 1925) the writer discussed the "energy- 
frequency spectrum" (hereinafter precisely defined) of irregular 
random disturbances extending over a long interval of time. In 
view of our lack of even statistical information regarding static or 
atmospheric disturbances the specification of the energy-frequency 
spectrum, denoted by R(co), was necessarily qualitative, and it was 
merely postulated that 

" R(co) is a continuous finite function of a> which converges to zero 
at infinity and is everywhere positive. It possesses no sharp maxima 
or minima and its variation with respect to co(co = 2irf), where it 
exists, is relatively slow." 

In a paper entitled "The Theory of the SchrotefTekt," J T. C. Fry 
deals with a similar problem, namely, the energy or "noise" absorbed 
in a vacuum tube from a stream of electrons with random time dis- 
tribution. His method of attack is widely different from that of the 
present paper. In a more recent paper on "The Analysis of Irregular 
Motions with Applications to the Energy-Frequency Spectrum of 
Static and of Telegraph Signals" (Phil. Mag., Jan., 1929), G. W. 
Kenrick, by making certain hypotheses regarding the wave-form of 
the elementary disturbances whose aggregate is supposed to represent 
static interference, and by applying probability analysis, arrives at 
explicit formulas for the "statistical" or "expected" value of i?(co) 
for a number of different cases. 

I 

In the present paper the statistical or "expected" energy-frequency 
spectrum R(u>) of random disturbances is investigated by a method 
which is believed to be somewhat more general and direct than that of 
Kenrick. 2 The results are applicable to the Schroteffekt, telegraph 

1 Jour. Franklin Inst., Feb., 1925. 

2 Kenrick's analysis is based on a formula derived originally by N. Wiener instead 
of proceeding directly from the Fourier integral. 
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signals and similar disturbances. The writer, however, concludes that 
their application to "static" or "atmospheric" disturbances is of 
questionable value owing partly to our lack of the necessary statistical 
information regarding such disturbances and also to the fact that they 
cannot be expected to have the "quasi-systematic" characteristics 
necessary to the application of probability theory. 

The energy-frequency spectrum of a disturbance, as the concept is 
here employed, will now be defined. Let a disturbance «f>(/) exist in 
the epoch ^ / ^ T and let 

F(*«) = C(w) + iS(a) 

= I <I)(t)e iu 'dt. (1) 

■A 

Then, as shown in my paper referred to above, 

- I \F(iu)\ 2 du = I Vdt. 
"Jo Jo 

The energy-frequency spectrum is defined by the equation 



1 
so that 



G(«) « Lim^l/^aOl 2 , (2) 



J 100 1 s*T 

G(co)dco = Urn ^ I &dt. 
o T ~* x L Jo 



(3) 



It is on this last equation that the physical application of the concept 
of the energy-frequency spectrum rests; namely, that it determines 
the mean square value of <£>(/), as the epoch T is made indefinitely great. 
Its principal application in electrotechnics depends upon the further 
fact that, if $(/) represents an electromotive force applied to a net- 
work of impedance Z(ia>), the mean square current T- absorbed by the 
network is given by 3 

P = Urn i f Pdt = f' ff") </ M . (3a) 

t-*«TJ q J Q |Z(*«)| 2 

We now suppose that the function or disturbance $(/) is composed 
of a number N of elementary disturbances; thus 

*(') = Zamcf> m Q - /,„), (4) 

3 A somewhat more involved formula gives the mean power absorbed. See my 
paper referred to in the first paragraph. 
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the wth elementary disturbance being supposed zero until t = t m . 
If we now write 

/«(*"«) = c m + is m = f <i> m {t)e iat dt, (5) 

Jo 

it is easy to show by the methods employed in my previous paper 
that 



I^MI 2 = Lfl» 2 !/™MI 2 



N-l N 



+ 21 £ a m a n (c m c n + s m s„) cos o>(/„ - /,„) (6) 

ro—1 n=m+l 

+ 2 £ 2- ama n (c m 5„ - s m Cn) sin o>(J„ - / m ). 

m=l n=m+l 



This is more compactly expressible as 

I^MI 2 = X>m 2 |/™M| 2 



1 



+ 2E E |a^n-/mM-/»(-^)e M '"-''" ) }RealParf (6a) 
m=l n=m+l 

Now, obviously, if the amplitudes a it • • •, a N and the wave form of 
the elementary functions fa, • •• , 4>n are specified, G(w) is uniquely 
defined and determined by the preceding formula. This, however, is 
not the case in the problem under consideration, where at best the 
functions are specified only statistically by probability considerations. 
Under such circumstances, when the problem is correctly set and 
sufficient statistical information is furnished for its solution, we 
introduce the idea of the statistical energy-frequency spectrum i?(w) 
defined as follows: 

The statistical energy-frequency spectrum R(u>) is equal to the weighted 
average of G(ui) for all possible values of G(u>), the weighting being in 
accordance with the probability of the occurrence of each particular 
possible value. 

For example, the statistical value of a function f(x lt x 2 , • • •, x n ), 
where the variables x u ■ ■ • , x n are defined only by probability con- 
siderations, is, in accordance with the foregoing definition, 

dx!pi(xi) ■ dx 2 p 2 (x 2 ) • ■ I dx n p n (x n )-f(x lt x 2 , ■ • •, x n ), 

00 J—x J — x 

where p m (x m )dx m is the probability that x m lies between x m and 
x m + dx m . 
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To apply the foregoing concept and definition of the statistical 
value of a function to the problem at hand it is necessary to suppose 
that the typical impulse f m (iu>) is a function of w and certain parameters 
Xi, X2, • • • , X„, and that these parameters are statistically specified by 
probability considerations. Thus we suppose that p m (\ m )d\ m is the 
probability that X,„ lies between X m and X,„ + d\ m . G(co) will then be 
a function of o> and Xi, X 2 , • • • , X n , the amplitudes a-\, • • • , a,y being 
regarded as parameters, when defined by probability functions. We 
then have, in accordance with the foregoing, 

R(u) = I </X,/>,(X) • I r/X 2 /> 2 (\,) 

%J — 00 \) — CO 

X ■■■ f d\ n pn(K)G(u t X„ X 2l • • • , X„). (7) 

J — DO 

II 

To apply the foregoing to the simplest possible case let us suppose 
that the elementary impulses are all identical; a A = a 2 = • • • a iV = 1, 
and that their distribution in time is purely random. With these 
assumptions it follows at once from (6) that 

*(*>) = - If(ico) |i + 2 • t |/(,- M ) 1 2 1 ~ gg « T t r _ >00 . (8) 

T IT 10- 1 

If /(*0) 7* 0, this has a singularity at u> = 0; however 

Lim -= I &dt = I R(<a)du 

^ ., (9) 

= v I (p-dt + v 2 J 0rf/ • 

Here v = A^/7" = mean frequency of occurrence of the elementary 
impulses. This formula is in entire agreement with Fry's results for 
the Schroteffekt (I.e.). 

To consider a somewhat more involved problem, we shall suppose 
that the durations of the individual impulses and their amplitudes are 
distributed at random. We further denote the probability that the 
duration of any impulse, selected at random, lies between X and 
X + </X by p(\)d\. Correspondingly, q(a)da denotes the probability 
that its amplitude lies between a and a + da. The durations and the 
amplitudes are then the statistically specified parameters. 

We now postulate that $(/) is an alternating series of impulses of 
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the same wave form; i.e. 

Hi) = L (- i) m Om4>m(t - tm), 

4>m(t) = 4>(t), < t < X m 

= * > Xm, 

2m = Xi + X 2 + • • • + X m _i, 

and we denote the mean frequency of occurrence, N/T, by v. 

Substitution in the preceding formulas and straight-forward opera- 
tions give 

12(c) = Z C 9 a*q(a)da • f |/(i«, X)| 2 £(X)rfX 
* Jo Jo 

plus the real part 

-\ (* aq(a)da\ • C f{iu,\)p(\)^d\ ■ f°° '/(- iu,\)p(\)d\ 

* |_ Jo _| Jo JO 

X Lim-^E £ (-l) n - m [ f%(X)^VxT l . (10) 

N—ras iV »n=l n=m-(-l L Jo J 

If we write 

r p{\)(**d\ = p(iu>) = P , (ii) 

we have by straightforward procedure 

1 N N r /"• "ln-m-l _ 1 

iis»£ J&.<- HI * X) H -****■ (12) 

whence 

R( a ) = 1 J" a*q(a)da j"° \f(io, X)| 2 />(X)JX 

7T |_ J * J I 1 + P(*w) J Real Part 

where 

/(*«, X) = f <p(t)e iat dt = c(co, X) + «(«. X), 

J7( w ) = P/Cito, X)/>(X)^MX, (14) 

V(«) = Cfi-ia, X)p(\)d\. 



STATISTICAL ENERGY-FREQUENCY SPECTRUM 379 

If, on the other hand, we suppose that the impulses, instead of 
systematically alternating in sign, are equally likely to be positive or 
negative, the double summation term of (9) vanishes and 

R(ai)=-f a-q(a)da- \ |/(ico y X) [ 2 /?(X)^X. (15) 

«" J-oo Jo 

This follows from the fact that the amplitude a is equally likely to 
be positive or negative. Consequently the integration with respect 
to da must be extended from - m to + w and, since by hypothesis 
q( — a) = q(a), it follows that 



f 



aq(a)da = 0. 



To apply the preceding formulas to actual calculations, it is necessary 
to know the function f(iu>, X) and in addition the probability functions 
involved. These latter may be supposed known from statistical data 
or calculable on theoretical assumptions. For example, if we assume 
that the times of incidence of the elementary disturbances are dis- 
tributed entirely at random, the application of well-known probability 
theory gives p(\) = ve~" x . 

A third case is of interest. Here, instead of postulating that the 
termination of one impulse coincides with the start of the next (i.e. 
tm+i — tm + X m ), we suppose that the times of incidence are entirely 
unrelated, and that the amplitudes are equally likely to be positive 
or negative. For this case the formula for i?(co) is formally identical 
with (15). 

Ill 

The foregoing analysis will now be applied to deriving what repre- 
sents more or less accurately the statistical energy-frequency spectrum 
of telegraph signals. To this end we shall suppose that the elementary 
disturbance may have any one of three possible values (all equally 
probable), characterized by durations Xi, X 2 , X3 and amplitudes 
o-i, o-z, 0-3- The corresponding spectra of the elementary disturbances 
are then determined by the equations, 

Miu) = f ' <t>(t)e iut dt, 
Jo 

/ 2 (*V) = \ <i>{t)e m dt, (16) 

Jo 

/s(*w) = f ' <t>(t)e«"dt. 
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The application of the preceding analysis to this case gives 

R( u ) = JL (ai f |/i(««) | 2 + a» a |/«(««) I 2 + fl3 2 |/ 3 M | 2 ) 

plus the rro/ par/ of 

X (ai/i(- *w) + a 2 / 2 (- ♦») + a 3 / 3 (- *«)) (17) 

XLim-f/z £ [*(** + **■ + **•)]-— l - 

JV— *■« xV ro = l n=rM+l 

It is to be understood that the raz/ par/ of the second term is alone to 
be retained. 
If we write 

1 I / N — 1*1— x^ -1 \ 

i £ £ [«**• + *■*■ + «"•)]-*-' =— ,(— " 7^-rr^) 

and 

}™*zz = xh x<0 

N 1 

= T *- L 

There is therefore an infinity at u> = 0, as we should expect. Its 
measure, however, is finite. 

The preceding is merely an example which admits of extension to 
more complicated types of signals, as will be obvious to the reader. 
For example, the probabilities of the elementary signals need not be 
the same and their number need not be restricted to three. 

IV 

In all the cases discussed above it will be observed that the dis- 
turbance is "quasi-systematic" in the sense that the elementary 
disturbances are all of the same wave-form differing only in duration 
and amplitude. Indeed, some such assumptions as these are essential 
to the application of the mathematical theory. In the case of atmos- 
pheric disturbances we have no reason to suppose any such quasi- 
systematic character exists. Furthermore, even if for the sake of 
argument, we suppose that the elementary disturbances, which make 
up static, have a common wave form at the point at which they 
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originate, they would vary widely in this respect after arriving at a 
common receiver. The writer is therefore of the opinion that the 
quotation from his previous paper appearing at the start of this article, 
represents all that can safely be said regarding the spectrum of static 
and that our present knowledge is insufficient to justify the application 
of probability analysis to the problem. All that we can say is that 
the part of R(u) which contributes to "static interference" is simply 



v 1 N 

Lim - . — £ a m 2 \f m (iu)\ 

AT-mi 7T iV 1 



a result deducible from (6) and in agreement with the conclusion of 
my original paper (I.e.). It is here supposed that the times of incidence 
are distributed at random. This formula, however, supplies no useful 
information in the absence of data regarding the wave forms and 
amplitudes of the individual disturbances. 



